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We study the relativistic quark-antiquark system embedded in magnetic field (MF). The Hamil- 
tonian containing confinement, one gluon exchange and spin-spin interaction is derived. We analyt- 
ically follow the evolution of the lowest meson states as a functions of MF strength. Calculating the 
one gluon exchange interaction energy (Voge) and spin-spin contribution (ass) we have observed, 
that these corrections remain finite at large MF, preventing the vanishing of the total p meson mass 
at some Bcrit, as previously thought. We display the p masses as functions of MF in comparison 
with recent lattice data. 
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I. INTRODUCTION 

During the last years we have witnessed an impres- 
sive progress of the fundamental physics in ultra-intense 
magnetic field (MF) reaching the strength up to eB ^ 
lO^^G ^ toJ [1]. Until recently magnetars [2| were the 
only physical objects, where such, or somewhat weaker 
MF could be realized. Now MF of the above strength 
and even stronger is within reach in peripheral heavy ion 
collisions at RHIC and LHC Q. High intensity lasers 
is another perspective tool to achieve MF beyond the 
Schwinger limit Q. On the theoretical side a striking 
progress has been achieved along several lines. It is be- 
yond our scope to discuss these works or even present a 
list of corresponding references. We mention only two 
lines of research which have a certain overlap with our 
work. The first one [a, Q is the behavior of the hydrogen 
atom and positronium in very strong MF. The second one 
[3| is the conjecture of the vacuum reconstruction due to 
vector meson condensation in large MF. The relation be- 
tween the above studies and our work will be clarified in 
what follows. 

Our goal is to study from the first principles the spec- 
trum of a meson composed of quark-antiquark embed- 
ded in MF. Use will be made of Fock-Feynman-Schwinger 
representation (see |8|] for review and references) of the 
quark Green's function with strong (QCD) interaction 
and MF included. An alternative approach could have 
been Bethe-Salpeter type formalism. However, for the 
confinement originating from the area law of the Wilson 



* andreichicov@ niail.rul 
t borisk@itep.ruJ 
-f orlovskii@itep.ru] 
§ ,sinionov@itep.ru| 



loop, the use of the gluon propagator is inadequate. Nu- 
merous attempts in this direction failed because of gauge 
dependence and the vector character of the gluon prop- 
agator, while confinement is scalar and gauge invariant. 
Therefore it is sensible to use the path integral technique 
for QCD-fQED Green's functions. This method based on 
the proper-time formalism allows to represent the quark- 
antiquark Green's function via the Hamiltonian (see Q 
for a new derivation), and it was used in [10] to construct 
explicit expressions for meson Hamiltonians without MF. 
In this way spectra of light-light, light-heavy and heavy- 
heavy mesons were computed with a good accuracy, us- 
ing the string tension cr, strong coupling constant as and 
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quark current masses as an input 

In what follows we expand this technique to incorpo- 
rate the effects of MF on mesons. The latter contains: 
1) direct infiuence of MF on quark and antiquark, and 2) 
the infiuence on gluonic fields, e.g., on the gluon propa- 
gator via qq loops and on the gluon field correlators de- 
termining the string tension cr. Since MF acts on charged 
objects, its influence on the gluonic degrees of freedom 
enters only via (iVc)^*', k — 1,2, ..., however, corrections 
of the second type can be important, as shown in [ij]. 3) 
As was shown recently in the framework of our method, 
MF also changes quark condensate (qq) and quark de- 
cay constants f^ etc., and in this way strongly influences 
chiral dynamics [141] . 

The important step in our relativistic formalism is 
the implementation of the pseudomomentum notion and 
center-of-mass (cm.) factorization in MF, suggested in 
the nonrelativistic case in [15| for neutral two particle 
systems. Recently the cm. factorization was proved 
for the neutral 3-body system in jig , the situation with 
charged 2-body system was clarified and an approxima- 
tion scheme was suggested in [T^]. 

The plan of the paper is the following. Section 2 con- 



tains a brief pedagogical reminder of how the two-body 
problem in MF is solved in quantum mechanics. The cen- 
tral point here is the integral of motion ( "pseudomomen- 
tum") which allows the separation of the center of mass. 
Here we also show how to diagonalize the spin-dependent 
interaction. In section 3 we formulate the path integral 
for quark- antiquark system with QCD-f QED interaction. 
Then from Green's function the relativistic Hamiltonian 
is obtained. Section 4 is devoted to the treatment of con- 
fining and color Coulomb terms. Here we also present the 
derivation of the eigenvalue equations for the relativistic 
Coulomb problem. In section 5 we discuss the spectrum 
of the system focusing on the regime of ultra-strong MF. 
Section 6 contains the discussion of the results, compar- 
ison with lattice calculations, drawing further perspec- 
tives and intersections of our results with those of other 
authors [l-Q]. 



II. PSEUDOMOMENTUM AND 
WAVEFUNCTION FACTORIZATION 

The total momentum of N mutually interacting par- 
ticles with translation invariant interaction is a constant 
of motion and the center of mass motions can be sepa- 
rated in Schroedinger equation. It was shown 15] that a 
system embedded in a constant MF also possesses a con- 
stant of motion - "pseudomomentum" . As a result for the 
case of zero total electric charge Q = the cm. motion 
can be removed from the total Hamiltonian. The sim- 
plest example is a two-particle system with equal masses 
TOi = m2 = m and electric charges ei — —e^ = e. We 
define 



R 



ri +r2 



J7 = ri ra, P = Pi+P2. 



(1) 



Straightforward calculation in the London gauge A 
i(B X r) yields 

One can verify that the following "pseudomomentum" 
operator F commutes with the Hamiltonian ([2]) 



F = P + -(Bx,7). 



(3) 



This immediately leads to the following factorization of 
the wave function (WF) 

*(R, T]) = ip{T]) exp jiPR - i-{B x j7)r| . (4) 

For the oscillator-type potential V{r]) the problem re- 
duces to a set of three oscillators, two of them are in a 
plane perpendicular to the magnetic field and their fre- 
quencies are degenerate, while the third one is connected 
solely with V{ri). 



Next we briefly elucidate the spin interaction in pres- 
ence of MF. The corresponding part of the Hamiltonian 
may be written as 



Hs = a,,/(cri(T2) - ^J,B{(Ti - 0-2), 



(5) 



where ei — —62 = e > and /x > 0. Diagonalization 
of Hs yields the following four eigenvalues e.g. for uu 
system, comprising both p and vr levels. 



^1% = o-hf, E^A - i"/!/ I 2l 



\ahf 



(6) 



where we assume that B is aligned along the positive z- 
axis and B = |B|. In a strong MF when jiB > Uhf spin- 

(s) 

spin interaction becomes unimportant and £"3 4 ~ ±2p,B. 

(s) 

For the lowest level El this corresponds to a configura- 
tion I H — ) when the spin of negatively charged particle 
is aligned antiparallel to B, and the spin of the positively 
charged one - parallel to B. This means that the spin 
(and isospin) are no more good quantum numbers and 
eigenvalues ^ correspond to the mixture of spin 1 and 
spin states. As a result the qq state will split into 4 



(s) 



r(-)^ 



states (two of them coinciding £[ = £"2 )■ Till now 
we treated a nonrelativistic system, to incorporate rela- 
tivistic effects we shall exploit the path integral form of 
relativistic Green's functions [3, Q • 



III. RELATIVISTIC qq GREEN'S FUNCTION 
AND EFFECTIVE HAMILTONIAN 



The derivation of the relativistic Hamiltonian of the qq 
system in MF consist of several steps. The first one is the 
4d relativistic path integral for the qq Green's function. 
The starting point is the Fock-Feynmann-Schwinger 
(world-line) representation of the quark Green's function 
3. The role of the "time" parameter along the path 
(si) of the z-th quark is played by the Fock-Schwinger 
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proper time Si,i — 1,2. Consider a quark with a charge et 
in a gluonic field A^ and the electromagnetic vector po- 
tential Ajf , corresponding to a constant magnetic field 
B. Then the quark propagator in the Euclidean space- 
time is 

S,{x, y)^{m, + d- igA - ie.^A^'^)-^ = (m, + D^'^)~^. 

(7) 
The path-integral representation for Si 3 is 



S,{x,y) = (m,-^«) / ds,{D^z),ye'''^^l:\x,y) ee 

= {m,~m)G,{x,y), (8) 



where 



Ki — m. Si + — dr, ' 



W^ 



dr. 



$^*' (a;, y) = PA^i^^ cxp ( ig / A^dz^'^ + ie, / A^f^dz^*) ) ^xp ( / dT,a^^{gF^^ + e^B^^,) 



(9) 



(10) 



Here F^i^ and J5^jy are correspondingly gluon and MF tensors, Pa,Pf are ordering operators, tr^i/ = jiilfj.1i' — lulfj.)- 
Eqs. ((71 ITU1) hold for the quark, i = 1, while for the antiquark one should reverse the signs of e^ and g. In explicit 
form one writes 



(tH crE 



'^^"-^'^'^ ^ 1 o-E crH ) ' '^'^'^^A'"' - 1 (tB 



rB 



(11) 



Next we consider qiq2 system born at the point x with the current jviix) = qi{x)Tiq2{x) and annihilated at the 
point y with the current jpa (y)- Here x and y denote the sets of initial and final coordinates of quark and antiquark. 
Using the nonabelian Stokes theorem and cluster expansion for the gluon field(see _11.] for reviews) and leaving the 
MF term intact, we can write 



OO /'OO 



"'0 



Gq,<iAx,y)= dsi ds2{D^z^^\y{D^z^^\y 



-K1-K2 



tT{fwAA))A>^ 



X exp (tei f A(f)dz(,i) - ie2 f ^'^d^'^ + ei / ' dnitrB) - 62 f ' dr2(crB)) , (12) 

\ Jy Jy Jo Jo / 



where 

f ^ Tiimi - Di)r2{m2 - D2), 
and Fi = 7^, F2 — ju for vector currents, F^ = 75 for pseudoscalar currents, while 

{WM))a = exp (-^ f dTT^,{l)dn)^„{2){F^,{l)Fx„{2)) +0{{FFF)) 



(13) 



(14) 



where dTr^i/ = dSf^^, - 



.(1) 



.(2) 



o'fj,JdTi — alJdT2, and ds^i/ is an area 
element of the minimal surface, which can be constructed 
using straight lines, connecting the points z^ (t) and 
Zj '(t) on the paths of gi and 52 at the same time t |8l.[l0|. 

Note, that operator T actually do not participate in field 
averaging procedure; as was shown in [l8|, the following 



replacement is valid: m~ D 



m - ip, p^ 



1 { dz^\ 

2 \dr J,^ 



As a result of the first step the qq Green's function is 
represented as a 4d path integral (including Euclidean 
time paths) and in addition also integrals over proper 
times Si, S2- In the second step one introduces monotonic 
Euclidean time tsiT) — X4 + -T, where T = \x4 — J/4I, so 
that 2:4 (r) = i£;(r) + Az4(r), where Az4(r) is fiuctuation 
of time trajectory around tEir). This new variable ts is 
an ordering parameter for trajectories z'^^^(tB), z'^^^{tE), 
and proper times transform into physical parameters ~ 



virtual q and q energies uji = o— , so that dsi — 



rduJi 



Combining for simplicity all fields into one Wilson loop 
W{A, A*^*^)), one can rewrite the Green's function in new 



I 

variables as 

Gqiq2{.X,y) 



T 



duji d(jj2 



LO 



3/2 3/2 



(i^3^(l)i^3^(2)) 



X e 



-Ki{bJx)-K2{u)2) 



{{fWF))^z„ (15) 



(see [9| for details of derivation). Here Ki{u>i), K2{u>2) 
are obtained from Ki in (j9]) by the same replacement 



dTi '- dtE ' 



Ki{uji) + K2{uj2) = I -^- + -^- 1 T~ 



dtf 



2uji 

wi / dz(i) ^ ' 
2 



dt 



E 



2W2 
tJ2 / dz(2) 



(16) 



The final step is the use of the Wilson loop dynamics 
to express all dynamics in terms of instantaneous inter- 
action. Indeed, the quadratic field correlator in (fT4|) is 
represented thr oug h two scalar functions D{z) and Di{z) 
(see, e.g., [ll|, [l^l for details), first of them is respon- 
sible for confinement, while the second one gives one 



gluon exchange (OGE) potential. So, for the case of zero in _D(j/, Zi — Z2), I?i(z^, Zi — Z2) one obtains a simple 

quark orbital momenta with the minimal surface, dis- instantaneous answer for spin-independent (SI) part of 

cussed above, integrating over relative time v = t\j — t'^ {Wa{A))A, 

I 



{W„{A))'X' =exp 



dt^ 



,(i)_^(2)i :!. 



a. 



3|z(i)-z(2)| 



(17) 



containing Vconti'i') = err and VbGE('') = ^^§p- Here a is the QCD string tension, cr = 0.18 GeV^ in our calculations. 
First we need to find the Hamiltonian Hg-^g^ of the system at if = i|' = t^. To this end we define the Euclidean 

Lagrangian Lg^g^. We write ^§:p- = 2uji-^^ = 2uJiZk, k — 1,2,3. Then all terms in the exponents in (IT^ . p^ and 
([T7)) can be represented as exp(— J dt^ Lg_^g^) and thus we arrive at the following representation: 



Gqiq2ix,y) = Z- 



T f f dujiduj2 

8^7 J (c^iC^2)3/2 




{D'z^'^D'z^'^)^ytr(^e-'t,2f'j 



(18) 



with the action 



'0 



/ ^t^ 2:(?(4")^— 4"4-V^H-^.^.=.,„, 



,{e)M)\ W1+W2 mf 



1-62- 

2u}o 



+ (t|zW-z(2)| 



as 



3|z(i 



5(2)1 



(19) 



Here A)^ is the fc-th component of the QED vector potential. The next step is the transition to the Minkowski 
metric. This is easy, since confinement is already expressed in terms of string tension. We have exp(— J L^dts) —>■ 
exp{i J L'^^dtui), tE — > JiAf, and 



Pk - 



M 

- f.j.^W _i_ p./l('') U - 



4 Pk 



(«)^(0 



^M 



(20) 



Explicit expression for Hamiltonian without spin-dependent terms is 



_ ^ (pW~e.A(zW))^+mf + c..re.^('^B ,,, ,,, 4 



as 



3|z(i)-z(2)|' 



(21) 



The qq Green's function ([TS]) takes the "heat-kernel" form, when going back to Euclidean time with Hamiltonian 



Gqiq2{x,y) — 3— 



diOl f diO 



8^./ a;?/^7 uj'^' 



tr(re"^«i 



(22) 



The cm. projection of the Green's function yields 

CO 00 

T f duji f duj2 



Gg^g^{x,y)d\x-v) = — 



8^./ ,,,3/2 



3/2 



Y,^l{Q)Mf)), 



-M„{uii,uj2)T 



^2 ri=0 



(23) 



r 



where ipn and M„ are eigenfunctions and eigenvalues of the effective energies w, are to be found from the mini- 
Hamiltonian Hg^g^. At large T the integral over U!i,uj2 mum of the total mass M„(wi,a;2), as it was suggested 
can be taken by the stationary point method, and hence in [1(|. To introduce the minimization procedure and to 



check its accuracy we shall begin by the calculation of 
the eigenvalues of one and two quarks in MF, and the 
energy of the ground state of a relativistic charge in the 
atom in the next section, reproducing the known exact 
results. 

We have the following equations defining Ui from the 
total mass M(a;i) 

7?^ = M„(c.i,c.,)^, ^^''f^'"^^ =0. (24) 

OUJi 

For a single quark in MF the first of the above equa- 
tions gives 

,,, , pl+ml + \eB\{2n+l)-eBa, tj 
A/„H^ 2^ +2- (2^) 

Then the minimization over uj yields the correct answer 

M„ = [pI + m\ + \eB\{2n + 1) - eSa^)!/^. (26) 

Now we turn to the case of q\q2 system and introduce 
the coordinates which are the generalization of ([T]): 



R 






aR' 



^ 



.(1) 



dr\' 



,(2) 



(27) 



(28) 



It is convenient to introduce the following two addi- 
tional parameters 



UJ\UJ2 



Wi - i02 



(29) 

cji + a;2 Wi + W2 

Let us consider the case of neutral meson, so that e\ = 
— e2 = e. Then the total Hamiltonian may be written as 

(30) 





H,,q, ^Hb+H^ + W, 


where 




2wi 


P + TT -Bx Rh 77 

a;2 2 V ^1 /. 


1 

^ 2UJ2 


— P-7r + -Bx R ry 

_Wi 2 V ^2 /. 



2w V 



6 C 

, TT B X R -I- s-B x 77 

2w V 2 2 ' 

1 



X P -B X 77 

2(wi+W2) V 2 ' 



(31) 



Equation (|3T|) is an obvious generalization of ([2]). The 
two other terms in (1301) read 



^ ^ ml -I- a;2 - ecriB _^ m| + a;f + ecr2B 



2wi 



VF = Fconf + VoGE + AVK = Cr?/ 



2a;2 

4 as (77) 
3 7? 



AVF, (33) 



and AW^ contains self-energy and spin-spin contribu- 
tions, which come from unaccounted spin-dependent 
terms of {Wa{A)). One can verify, that the "pseudomo- 
mentum" operator in (I3|), introduced in section II, com- 
mutes with Hb and hence we can again separate the cm. 
motion according to the ansatz (|4]): 



i/B*(R,r7) 



exp < iPR— i-(B X 



77)r} 



HB^piri). 

(34) 

Then the problem reduces to the eigenvalue problem for 
'^{rf) with the Hamiltonian Hb having the following form: 



1 



Hb = — \ -i^ 



2w 



d 



dr] 



s-B X 77 



2(cji + W2) 



(P - eB X 77)^ (35) 



For P X B = the system has a rotational symmetry 
and the cm. is freely moving along the 2-axis. Here 
we shall consider a state with zero orbital momentum 



(Lr,). 



^ ^ idT] 



0. As a result Hb is replaced by 



a purely internal space operator 



i^o 



1 

2(^ 






-(Bx77) 



(36) 



To test our method we put W = ^ and arrive at the 
equation 



(ifo + ha)ip = Af(iui,W2)'^- 



(37) 



Consequent minimization of Miuji^txiq) in a;i,W2, as in 
(l26t , yields the expected answer for the two independent 
quarks. 



M = Jm{ -H eB(2ni + 1) - ecriB- 



i^ -t- eB(2n2 -t- 1) -t- e(T2B. (38) 

We turn now to the particular case of charged two- 
body system in MF, ei = 62 = e and also ttii = ?7i2, 
when exact factorization of R and 77 can be done. In this 
case, for wi = a;2 = w and P x B = 0, the Hamiltonian 
has the following form [9| 



2m2 + 2ij? - e{(Ti + cr2)B a ( rf 

^ ^2 



2uj 



7 



l^OGE + Vss + AAfsE. (39) 



IV. TREATING CONFINEMENT AND GLUON 

EXCHANGE TERMS. THE ABSENCE OF THE 

MAGNETIC QCD COLLAPSE 

From (P5| . (|55| it is clear, that inclusion of T4onf and 
Voce in ^0 + W leads to a differential equation in vari- 
ables i^^jTjz, which can be solved numerically. However, 



in order to obtain a clear physical picture, we shall rep- 
resent Vconf in a quadratic form. This will allow to get 
an exact analytic solution in terms of oscillator functions 
with eigenvalue accuracy of the order of 5%. The OGE 
contribution will be estimated as an average ((^JVogeI'/'), 
thus yielding an upper limit for the total mass. 
For Vconi we choose the form 



K. 



V,, 



a I rj 



1 



+ 7 



(40) 



Here 7 is a positive variational parameter; minimizing 
Vconi w.r.t. 7, one returns to Vconf- We shall determine 



M(a;iW2,7) corresponding to 
additional condition 

i9M(wi,W2,7) 



Konf, and to define 7 an 



97 



= 



(41) 



Minimizing in wi in the limit 7712 3> ttii (the hydrogen 
atom), one obtains 



M = mii/l - a^ 



TO2, 



(44) 



which coincides with the known eigenvalue of the Dirac 
equation. 

In our ((7192) case one can calculate the expectation 
value of VoGE — — | " with the asymptotic freedom 
and IR saturation behaviour in p-space (see [2^ for a 
derivation and a short review) 



asiq) 



An 



/3oln(4^) 



(45) 



will be added to dM]) ■ As a result M{Jf^ , 0;^°^ , 70) will be 
the final answer for the mass of the system, neglecting the 
AW contribution. The difference of the exact numerical 
solution from that obtained with the genuine potential 
Kont does not exceed 5%. The solution of the equation 
(-ffo + Kont)^' — A/(a;i,a;2,7o)</' for the ground state is 



V'(^) 



1 



71-3/2 j,2^ J, j^ 



: exp 



2r\ 



2rl 



(42) 



where 



rA. 



eB \^ ^ -ye'-'B'' J 



-1/4 



ra 



(— ) 



1/4 



As 



we shall see below, for the lowest mass eigenvalue with 
eB ^ cr, one has r_L ~ W^, tq ~ -^ and the (9192) 
system acquires the form of an elongated ellipsoid. Sim- 
ilar quasi-one-dimensional picture was observed before 
for the hydrogen-like atoms in strong MF [3, Q . In such 
geometrical configuration Voge manifests itself in a pe- 
culiar way, again similar to what happens in hydrogen, 
or positronium atoms, and as was shown in jl3| in QCD 
the outcome is also similar to the case of QED, with the 
screening of the diverging effects. 

We turn now to the OGE term to treat it in our for- 
malism. As a starting point we present another check of 
our approach, namely we shall obtain the ground state 
energy of two relativistic particles with opposite charges 
without MF interacting via the Coulomb potential. The 
corresponding Hamiltonian reads H = Hq 
H<f) = M(l), and for eB = we have 



i/,-^,then 



M ■ 



2uji 



ml + UJ2 

2UJ2 



(43) 



where /3o = ^-iV e — f"-/: ^^b is proportional to ^/a, 
Mb w 1 GeV ^. With the wavefunction ^^ the aver- 
age value of Fqge takes form 



AM, 



OGE 



VbGE(<7)V'^(q) 



dh 



4 
37r7o 



' (27r)3 
as{q)dqe ^ / 



q^{rl-rl) 



(46) 



where ip'^{q) is the Fourier transform of squared wave 
function ilP'{r}) and I{a^) = j_^ dxe^"' ^ . Estimating 
the integral in (|46)) . for eB ^ a, i.e. for tq ^ r± one 
obtains for massless quarks 



AM, 



OGE 



3ro/3o 



In In ■ 



/ct In In ■ 



eB 



(47) 



With eB increasing the upper bound for the qq mass 
is boundlessly decreasing. The exact eigenvalue should 
lie even lower. 

This situation is similar to the hydrogen atom case, 
where AMqouI diverges as — In eB, and in this case e~^e~ 
loop contribution to thephoton line stabilizes the result 
(the "screening effect" [y|, H)- In our case the qq loop 
contribution to the OGE term can by written in a similar 
way, adding to the gluon loop also the Lowest Landau 
Level (LLL) of the qq in the MF, 



VbGE(Q) - - 



IBttq:, 



(0) 



47r 3 ^ 



Mo 



ai°'nj|e„B| 



exp I 



W7b\ 



T 



m 



(48) 



where 




10* 10° 10° 
|eB|, GeV' 



FIG. 1. One gluon exchange correction to the meson mass in 
GeV as a function of magnetic field with (solid line) (|49|l and 
without (broken line) H46p account of quark loops contribu- 
tions. 



1 m/ \ ln( Vl + '2 + \/^) -1 i^ 1 1 • 1 

where 1 iz) — ^ , ' + 1. Calculating now the 

y/z{z + l) 



average value of 



AM, 



OGE 



(Vc 



oge; 



(49) 



one obtains saturation of AMqge at large eB, as shown 
in Fig. 1, eliminating in this way the possible "Color 
Coulomb catastrophe", discussed in the first version of 
this paper |23j . 



V. MESON MASSES IN MAGNETIC FIELD 

Our next task is to calculate analytically the mass 
Mn{uJi,uj2,7) of a (5192) meson. We have to solve the 
equation 

(Ho + H, + W)^niv) - Mn{iOi,i02n)-^n{v), (50) 

where Ho,H„,W are given in ^, ([33]), ([36]), the total 
Hamiltonian for charged meson is given in p9|) . 

The resulting mass for neutral meson without spin- 
dependent contribution from AW is 



+ ■ 



M"„(cji,i:j2,7) =£: 
m,2 + UJ2 + eB(T2 

2U}2 



-AM, 



OGE- 



M„(wi,W2,7)^ 



■ml+iul — eB<Ti 
eB<Ti eB<T2 



2^1 



2^2 



(51) 



2a; 



'e2B2 



4(7(1; 



-(2n_L + l) 




?. (52) 



AMqge is given by (05]) and P^ . So, for fixed n we 
have four states for different quark spin orientations, | -|- 

-I-), I -I — ), I — h) and I ), where +/— are up/down 

directions of individual quark spins, with corresponding 
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(53) 
(54) 
(55) 
(56) 



The spin-dependent part AW contains self-energy Vy 
and spin-spin Vss contributions. As was shown in |24| 
the mass correction, corresponding to Vse, is given by 



AM: 



SE 
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A'KLJl 



1 + 7/ \\2eB 
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47ra;2 
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■n 



\x 2eB 



(57) 



where ri{t) — t J^ z^Ki{tz)e''''dz and A ^ 1 GeV^^ is 
vacuum correlation lengths. 

Let us introduce now the spin-spin interaction. It has 
nondiagonal structure 



Vss 



STTtts 

9a;ia;2 



,5(3) (: 



r)(Ti(T2 = a55cricr2. 



(58) 



so we should diagonalize the Hamiltonian with respect 
to spin variables. This results in new four states, two of 
them are mixture of |-| — ) and | — h) states, corresponding 
to 7r° and p" with zero spin projection Sz = 0, the other 

two states | + +) and | ) correspond to p° states with 

Sz — ^ and Sz — —1 (we consider ground state n = 0). 
We note that both Vss and AAIse are to be considered 
as corrections and contain ojl , ^2 ? obtained from min- 
imization of the remaining part of the meson mass. Note 
also, that masses M^ and M~+ are symmetric with 
respect to uji -H- uj2 for equal quark masses, so in this 
case we have in fact only two variables in minimization 
procedure w and 7. 

The masses of first two states are: 



^1,2 = ^(Mn 



M22) ± 



'M: 



22 



Ml] 



4ai2a2i, 
(59) 



where 
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(60) 



_ (o)_, , and {ass) is the aver- 



ai2 = 021 = \ass/|, ,(0) 

■^1 - - , 

aging with the wave function (j42|) (see [24] for derivation). 
The parameters a;± and Wip are obtained by minimizing 
of corresponding diagonal eigenvalues Mq and Mq~^, 
the parameter 70 (see (35)) for the ground state is de- 
fined from the condition 



dMo{uji,uj2,^) 



97 



97 



0. 



(61) 



It is easy to see, that at large eB masses -Ei^2 tend to 
diagonal values 

£■1(6^^ 00) ^Af22, E2{eB -i' 00) -> Mil. (62) 

The remaining two states have masses 

E3 = M+++AMsE+{ass), Ei = Mr, 



,(0) ,,(0) .,0 



-AMsE+{ass), 
(63) 

taken in point {uj''^\uj2'' ,'y'^) in accordance with mini- 
mization conditions (P^ and (pTj) . 

It should be noted, that actually we have eight states 
instead of four, since qq systems with different quark 
charges behave differently in MF, as we see from our 
Hamiltonians. Isospin is not conserved now and each 
neutral state splits into two states with different quark 
content uu and dd. 

Let us consider also the particular case of charged me- 
son with Hamiltonian (159]) in state with Sz = 1 {\ + +)- 
state, corresponding to p^). The eigenvalue, correspond- 
ing to this state, is given by the following expression 



Mn{uJ,-f) 




AM, 



OGE 



AM: 



SE 



{ass)- (64) 



Among considered states, the mass of charged meson 
ground state (p^ with Sz = 1) and -©2, corresponding 
to 7r°, tend to finite value at large MF due to cancella- 
tion of linearly growing terms in en^,n^ and in iJ^, while 
other masses grow with eB. This is true, provided that 
the spin-spin contribution {ass) remains finite at large 
MF. However, it contains the factor V'^(O) ^ eB, which 
leads to unbounded decrease of £'2- As was shown in [2J|, 
this situation is not physical, the total mass eigenvalues 
should be positive, and the reason of this decrease is the 
unlawful use of the perturbation theory for the poten- 
tial c^'^^^(r). One should replace ass by a smeared out 
version, e.g.. 



5(3)(i.)^^(3)(i.) 
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,-rVA^ 



A - 1 GeV" 



Using the wave function (|42|) . one obtains for {ass) 
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(66) 



The smearing length A on the lattice corresponds to the 
lattice unit a (A ~ a), in physical situation the relativistic 
smearing is connected with the gluelamp mass parame- 
ters in D{z) and Di{z), see [1^ for details. 



-A- M(p° j, s^=0 


' 1 ' 

A Hidaka, p° 


1 ' 


1 




1 


. -v- M(p°;, s^=0 


T Hidaka, p* 




-A--' 




_ 


-A- M(p°J, s^=1 












_ -v-M(p°;,s^=1 
M(p*), s^=1 


.-A'''' 
.A''' 




,--V' 




--^ 


.A-' 
.A-'' 


_...V'" 


-v" 


^ 





-I 


A'' .-V 

a .V' ^y^ 


--^^A 


A 


_^ 


____A__ 


t 




▼ 

1 


1 




▼ 




1 




■ 


1 1 1 



(65) 



eB, GeV' 



FIG. 2. The masses of the systems in GeV as a functions of 
eB. See the text for explanations. 

In Fig. 2 we plot the masses of some selected systems 
as a functions of eB (e is the p+ charge, not the charge of 
individual quarks). Calculations were performed accord- 
ing to dSn]), (|55)) . ([M)) and the minimization procedure. 
The dashed curves correspond to the p° state with Sz = 
(eigenvalue Ei), the solid-symbol lines describe p" state 
with Sz = 1, the lower solid curve refers to the state 
of charged meson p"*" with Sz = 1- The black triangles 
are from lattice calculations [25 • One can see that the 
masses of first states are increasing, while the last one 
tends to finite limit in accordance with discussion above 
(note, that the results plotted in Fig. 2 were obtained for 
massless quarks). 



VI. DISCUSSION AND CONCLUSIONS 

In our treatment of relativistic quark-aniquark system 
embedded in MF we relied on pseudomomentum factor- 
ization of the wave function and relativistic Hamiltonian 
technique. The Hamiltonian for mesons in MF, contain- 
ing confinement, one gluon exchange and spin interaction 
was derived. Using a suitable approximation for confin- 
ing force we were able to calculate analytically meson 
masses as functions of the MF. In this paper to simplify 



things we started with p" meson states at S = taking 
7i in place of Fi and T2 in p^ . In this way we essentiahy 
left aside the complicated problem of chiral dynamics and 
pseudo-Goldstone spectrum. In this oversimplified pic- 
ture the lowest neutral state with Sz = is a mixture of 
the p^ and tt^, as can be seen from its spin and isospin 
structure. Indeed, uu or dd system under consideration 
is a mixture of isospin 7 = and 1=1 states, and at 
large MF it has a spin structure \u 'f,u l), which is a 
mixture oi S = and 5* = 1 states. We have calculated 
mass of the higher state of this mixture, which we call 
p^{eB), while the lower state, associated with n^i^eB), 
can be subject to chiral corrections. These states have 
negative corrections from one-gluon exchange and spin- 
spin interactions. As was shown, these corrections (and 
the total mass) stay finite at large B, preventing the so 



called "magnetic collapse in QCD", discussed earlier in 

As shown in Fig. 2, our analytical results are in agree- 
ment with lattice calculations [25| both for p° and p+ 
states. 

Another system which can be treated using the same 
technique is the neutral 3-body system, like neutron. 
The results might be important for the neutron stars 
physics. 
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